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The  Analysis  of  Lock  Foundation  Design 


Introduction 


The  design  of  structures,  such  as  locks,  footings  and  culverts 
resting  on  a  semi-infinite  body,  is  dependent  on  the  analysis  of  the 
contact  pressure  between  the  slab  and  the  foundation.  For  footings  and 
culverts  the  assumption  normally  made  is  that  the  contact  pressure  is 
uniform  over  the  area  of  the  slab.  In  the  case  of  locks  and  similar 
strip  foundations  it  is  assumed  that  the  deflection  is  proportional  to 
the  pressure.  This  assumption  gives  rise  to  the  hypothetical  concept 
of  a  modulus  of  the  foundation.  By  means  of  this  simplifying  assumption 
the  contact  pressure  can  then  be  determined  from  a  solution  of  the  differ¬ 
ential  equation  (1)  subject  to  the  particular  boundry  conditions. 


-  ~  Q) 


a  7^  v 

q^) is  the  load  on  the  slab,  and  D  the  flexural  rigidity.  The  modulus 
of  the  foundation U  is  defined  by  (2)  as 


in  which  p  and  uo  are  the  contact  pressure  and  the  deflection  respectively 
for  any  given  point  of  the  slab. 

The  assumption  of  a  constant  modulus  of  the  foundation  is,  in  general 
incompatible  with  the  mechanical  behavior  of  solids  and  soils.  To  have 
reality  agree  with  theory  it  would  be  necessary  that  the  supporting  media 


« 


be  made  of  identical  but  individual  springs  each  of  which  could  function 
independently  of  the  other.  In  other  words  the  hypothesis  implies  that 
the  deflection  is  independent  of  the  loaded  area  and  that  a  load  can  be 
the  cause  only  of  the  deflection  underneath  the  loaded  area.  Computations 
based  on  the  use  of  equation  (1)  are  at,  the  most  crude  estimations. 
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The  theory  of  the  bending  of  a  bar  on  an  elastic  foundation  in 
which  the  hypothesis  of  the  modulus  of  the  foundation  is  adopted  dates 
back  to  1367..  In  1338  ^  the  method  was  used  in  the  design  of  railways. 
Mathematical  expressions  developed  from  (1)  for  the  case  of  concen¬ 

trated  loads  are  applied  to  the  design  of  highway  and  runway  slabs. 

For  such  applications  of  equation  (1)  it  has  been  shown  that  the 
computed  behavior  of  the  slab  checks,  with  fair  agreement,  field 

measurements.  The  U  value  for  such  designs  is  computed  from  field 

3  4 

tests  or  from  theoretical  results 

For  foundation  slabs  subject  to  uniform  or  varying  load,  over¬ 
turning  moments  or  concentrated  loads  at  the  edges  of  the  strip  such  as 
occur  in  locks,  +here  appears  to  be  no  obvious  justification  for  the 
vigorous  application  of  a  theoretical  analysis  based  on  the  hypothetical 

Ik  value.  Nevertheless  in  the  literature  there  can  be  found  many 
5 

examples  of  this  method  of  analysis.  The  modulus  of  the  foundation 
is  generally  determined  from  empirical  expressions  whose  limitations 
or  restrictions  are  meagrely  described. 

The  factors  which  affect  the  variation  of  the  contact  pressure 
are  the  size,  shape,  elastic  properties  and  flexural  rigidity  of  the 
slab,  the  elastic  properties  of  the  soil,  and  the  character  of  the 
superimposed  load.  Whenever  equation  (1)  is  to  be  adopted  for  a  pre¬ 
liminary  analysis  of  the  contact  pressure  beneath  a  slab,  a  good, 
approximation  for-  a  Ik  value  could  be  made  from  the  proper  selection  of 
pressure  and  deflection  values  as  determined  by  the  theory  of  elasticity. 


1 


Refer  to  the  list  of  references  in  the  bibliography 


. 

’ 


, 


:  ■-!  j 


.a\-  3  T  *  - 


' 


■ 


' 


■ 


3 


An  average  W  value  secured  between  the  limits  of  a  perfectly  flexible 

and  absolutely  rigid  slab  would  be  commendable.  Considering  a  few 

examples  we  have,  for' the  case  of  the  perfectly  flexible  circular  slab 

fen 

subject  to  a  uniform  load,  that  the  modulus  of  the  foundation  is  given  by 


U. 
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where  L(Jo  is  the  complete  elliptic  integral  of  the  second  kind,  a  ,  v 
and  E_  are  the  radius,  ^oisson’s  ratio,  and  the  modulus  of  elasticity 
respectively  for  the  slab.  It  can  be  seen  that  tv  depends  on  the  size 
of  the  slab  and  varies  from  a  value  of 


w 


-  -  -  (fr) 


at  the  centre  to  a  value  of 

N 


L 


at  the  edge.  On  the  basis  of  an  average  slab  deflection  the  value  is 

JL  'b  L  -  - —  (£>) 


wi 
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In  a  similar  manner  for  an  absolutely  rigid  circular  plate  the 

<ob 


modulus  of  foundation  s  given  as 

^  s  £ 


^  J  q1-  r1" 


which  gives  a  minimum  value  at  the  centre  of 


\H 


E. 


---W 


- - (8) 


to  infinity  at  the  edge 
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For  a  flexible  square  plate  the  W  value  computed  from  an  average 

deflection  is 

U..  *  OpBkb  It  -  -  -  W 

and  for  rectangular  plates 

~  —  —  0°) 
m  XJ 
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where  a  is  half  the  side  of  the  square  plate  and  f\  the  area  of  the 
rectangle.  Values  of  m  for  various  ratios  of  the  sides  of  the 
rectangle  are  given  in  reference  6a  page  333.  In  the  foregoing  computations 
it  has  been  assumed  that  noth  the  slab  and  "he  foundation  are  elastic 
isotropic  materials. 

Object 

In  the  first  part  of  the  work  to  follow,  consideration  will  be  given  to 
the  two-dimensional  problem  of  an  elastic  strip  resting  on  an  elastic 
isotropic  foundation.  The  solutions  will  be  carried  cut  for  concentrated 

loads  and  overturning  moments  at  the  edges  of  th©  strip.  The  me’ hod  of 

7 

analysis  will  be  similar  to  that  given  by  Borowicka  for  the  case  of  the 
uniformly  loaded  slab.  The  solution  takes  into  consideration  the  physical 
properties  of  the  slab  and  the  foundation.  No  simplifying  assumption 
that  the  deflection  is  proportional  to  the  pressure  is  made. 

The  second  part  will  include  a  discussion  and  criticism  of  the 
validity  of  Borowicka’s  method  of  analysis  for  the  cases  of  symmetrical 
but  discontinuous  loadings. 

In  the  third  part,  consideration  will  be  given  to  the  inverse 
problems.  That  is,  the  equation  for  the  foundation  reaction  is  assumed 
and  the  deflection  is  then  computed  from  Boussinesq's  equation  for  a  serai- 
infinite  loaded  mass.  Substitution  of  this  computed  deflection  into  the 
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5. 


slab  equation  (1)  yields  the  expression  necessary  to  determine  the 
superimposed  load  on  the  slab.  Examples  showing  the  determination  of 
the  contact  pressure  from  assumed  deflection  curves  will  also  be  given. 

Assumptions  and  Notation 

The  following  assumptions  are  made  with  regard  the  determination 
of  the  contact  pressure  and  deflection  for  the  slab  and  foundation; 

(a)  the  strip  or  slab  and  the  supporting  foundation  are  homogenous 
elastic  and  isotropic  materials. 

(b)  TIooke*s  Law  and  the  principle  of  superposition  are  valid. 

(c)  that  no  tangential  forces  exist  between  the  slab  and  the 
foundation. 

The  notation  used  is  as  follows; 


E  the  modulus  of  elasticity  of  the  slab 

s 

the  modulus  of  elasticity  of  the  foundation  assumed 
constant  for  the  full  depth  of  the  soil 
m_  Poisson’s  constant  for  the  slab 
%  Poisson’s  constant  for  the  foundation 
a  half  the  width  of  the  slab 
h  the  uniform  thickness  of  the  slab 
q  the  uniform  load  on  the  slab 
W  the  concentrated  line  load  at  the  edge 
M  the  overturning  moment  at  the  edge 

w  the  deflection  of  the  slab  or  foundation  positive  downwards 
x,  y  Gartesion  co-ordinates 


^  variables  in  the  x,  y  plane 

r  distance  from  the  load  to  point  at  which  deflection  is  computed 


. 


. 
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D 


the  flexural  rigidity 
of  the  slab 


N 


-l) 

=  mh4  Lb 


2_D_ 
tv\  a* 


a  constant  from 
Boussinesq’s  equation 
a  constant 


A0>  A2» 

B0,  ^2* 

c0,  cx, 

<<-  and 


A4 

b4 


<? 


) 

) 

)  coefficients 

) 

) 


-  constants 


Theory 

In  the  following  method  of  analysis  it  is  assumed  +hat  the  contact 
pressure  between  the  slab  and  the  foundation  can  be  represented  by  a 
Maclaurin  series.  With  this  assumption,  the  deflection  of  the  slab  is 
then  determined  from  the  slab  equation: 

1") 


d' 


oO 


and  the  deflection  of  the  foundation  can  bo  found  from  Boussinesq’s 
.  <o  b. 

expression 


00  =  X  -  l 


r 


- (a) 


For  the  deflections  as  determined  by  (11)  and  (12)  to  be  compatible  the 
coefficients  of  all  like  powers  of  x  in  the  two  equations  have  to  be 
equal.  In  this  manner  an  infinite  set  of  equations  involving  the  co¬ 
efficients  of  the  assumed  power  series  for  the  contact  pressure  are 
obtained.  The  contact  pressure  can  then  be  determined  to  any  desired 
degree  of  accuracy  from  the  solution  of  a  finite  number  of  the  infinite 


set  of  equations 
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Case  1 
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Concentrated  line  loads  at  the  Edges  of  the  Slab. 


F  l<q  1- 

For  Case  1,  shown  in  figure  1,  we  have  for  the  boundry  conditions 
of  the  upper  surface  of  the  foundation  that 

7^  y  K  -  o  uj  av\»^  Vq\u«s. 

^  =  jp(?0 


7s  <  Q 


The  boundry  conditions  for  the  slab  are 

=  q 

*  =  O 

The  slab  equation  for  case  1  is 

cT^o  =  - 

d 


a; 


-  0 


'  ~  (^b) 


uo  =  vjO  o  ^  ^  —  O 

4^  df* 


D 


-  -  -  (\5) 


The  solution  of  (15)  is  given  as 

VjJ  2  \jJ  ^  'A)  ^ 


where  w  is  the  solution  to  the  homogenous  differential  equation 

w 


d  ^  ^0 

a 


-  o 


-  -  Q\  fe) 


and  Wp  is  a  particular  integral  in  the  non-horaogenous  differential 

equation  (15). 

We  have  from  (16)  tha+  the  complementary  solution  is 

+  cqq)  *  cq*f  + 

The  particular  integral  is  assumed  to  be  a  Maclaurin  series  of 
even  powers  with  unknown  coefficients,  •‘•■hat  is 

OJ. 


—  In) 


=  £  ^  id 


A 


-  - 
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Likewise  it  is  assumed  that  the  contact  pressure  |^) is  represented 
by  +he  power  series 

rff"  -  — w 

0 

Substituting  (19)  in  the  slab  equation  (15)  and  inte^rqtincj 
we  have 


Comparing  the  coefficients  between  (13)  and  (20)  we  find 


X  q 

D 


b 


2.0 


Xr\\ 


The  deflection  of  the  slab  can  then  be  written  as 


i  *  4  I  r 

oj  =  +  Cjjc  \  -f  C*  (jjO  t  ''b  j jX  )  —  J-  7  \^xy\  uX. 

W  °  (*vu+)\  W 

From  the  boondry  conditions  of  the  slab  we  determine  that  the 


coeff ic ient 


Co  :  Oo 

C_t  r  C-3  =0 

c,=  V  .4 

and  the  deflection  for  the  slab  becomes 


10  =  U)0 


Z  bin  .  PL? 

X'O  yt";  0  (2.Y\t^Xntl)  ^ ^ 


- 
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V\5l> 


lY) 


ft 


(xv\t4) 
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X’A+d 


-  -  AO 


Considering  now  the  deflection  of  the  upper  surface  of  the 
foundation  which  must  conform  to  the  deflection  of  the  3lah  for  a  contact 
pressure  of  p(x),  we  have  that  the  deflection  w  at  any  point  x  due  to 


load  b's'y  is 


U 


"4 


1  i. 

1  ^  Y 


4''.  •  h\.  d\ 

4-S-y, 


>-0  5  2. 


■  V-;. 
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In  equation  (22)  let 

N  = 


and  replace  r 
variable  ^ 


we  have, 


<  v> 


Considering  first  the 


10. 


Since  duo  becomes  large  when  ^  -»=*=>  we  consider  the  difference  in  the 
deflection  between  any  point  x  and  the  origin,  thus  we  obtain 


-  c\vxj  5  |||  i  0 ^  +  \  t  '  -V  lc~  j  (  'Aj  “  1  l  | 


In  the  limit  as 

Itm 


^  the  first  form  reduces  to 

— n 

r  I 


\  0  cj 


\  i- 


»Prf 


\  oq  \ 


O 


i- 


Inte  -rating  with  respect  to  the  variable  *  we  then  have 


-  uo  ■=  X_ 

1\H 


|  \oc^Vd  ~ 

-ok  L 


ci - (*1) 


i  i  |  ^  +  i 


pjj) 


~1 


To  evaluate  the  first  two  integrals  in  (24)  we  define 

nY’  I  ~  P^-  -t  ( V“‘- 


-  -H 


(2n  \\)  Q 


xn 


(V 


and  F(}s)  =  O 


11. 


Integrating  lv  by  parts  we  have  that 

V  =  \  lo  j  l  *  - \) 


~  a 


'd  > 


-  a 


PiSiil 

Y*  ' 


-  *) 


Also  for  we  have 

1 4  1  ^  fY  10 


r^)d\ 

$-*■  ' 


1,  • 


jF0\) 

J  > 

K?)  loci  (.a -a)  -  [  F^)  d*j 

l  V"- 


Adding  \k  and  lv  we  find  that 

I\  tlj.-  FLM  Icq  U  - *)  -  r (r <0  ioq  u  +  ~  vvlL_  J  <\ . 

^  y.  ~  ^ 


»  ^ 

\y\tl  |  .  X^+f  m^1 

l  ocj  (,*-*)  d-  a _ 

,  m 


>  t  _  •  _  \  J-V^T1  \ntt  xn+i  in'” 

1,^1^  ■=  p7  E>i/\  \  a  -  /_  \  o c|  (jsv -  t*-)  -t  a  4-  >  locjU-t*) 

\  a2" 


ft  a  0 


Expanding  the  last  term  in  (25)  we  3ee  that 


(cV,'+‘  x*" 

-  ±  1  --  ^ 
^  j 

•  a  V 


-a‘VS't'  xvY> 

:i  -  >  . 

V" 


2n  an-r 
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^=o 


--Y) 
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which  upon  substitution  in  (25)  gives 
I»  +  I*  =  Bin 


1"+'  2.Y>-V\ 

a  -  * 


V\  *  0 


(j.  A  t  \)  h 

+  aXA  *'+>*"*)  \°^UtTt)  -  * 


v's  0 


\  0  t  C\  -  >) 


y  i  /  +'  Xn-V  "V  \ 

a  -  (r  *) 


Zn- r t \ 


From  the  last  Integral  ij  in  (24)  we  obtain 

l*  »  ■ 


!Z 


Etfc  n  

X.VV  t  ^ 


a  i  a 
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Adding  and  cancelling  we  have 

<p 

Vis)  o  "  '-O  -  Z. 

^  A* 9  X 


&xa.  ^  \  Uq  l  \-  JL\  + 


ll'-*)  ♦  m»  +  * 


-'!c^ 


i  ti 


-bL\  I  10 


X.V* 


Xrt  —  X  v"V  \ 


V=» 


Expressing  the  logarithmic  terms  by  their  series  expansion  we  write 


VO  5 


X. 


B 


x  v\ 


<X. 


i\~  o 


\aW 


<t> 

£ 


Bn 


SL 

r  U 


ti  mr 

VQ  1 


SA+V 


r=  i 


H 

r 


i-  PM 

v  q . 


Zntl  «° 


r=  i 


-L  (£ 

r  , 


r-\ 
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Adding  even  and  odd  powers  the  expression  for  the  deflection  of  the 


foundation  reduces  to 


A 


\ 

lf\-  2.V  + 1 


Lettering  Y-  f-l  and  then  rewriting  /  instead  of  r  ?;e  have 


o  -  xo 


vr 


(xr  f  X) 


JL 

a, 


xv\  -t  xf  i  2.  . 


V~3 


_  &  xn 

\*-r  ^xKx<ri  *0 


(^X^  -xr-  +1 


14. 


The  first  summation  expression  in  (26)  can  he  rewritten  in  terms  of 
one  variable  index  by  summing  along  diagonals  first  from  to 

and  then  from  to  as  shown  in  figure  3. 


F  \Ca 

The  final  expression  for  the  deflection  of  the-  upper  surface  of  the 
foundation  due  to  the  load  is 

=  J:  z.  -B—  -  It) 


r=t>  6 


(Tir-vaO^vx^) 


Expression  (26)  for  the  deflection  of  the  foundation  agrees  with 

Borowicka’s  t 
from  (27)  as 


Borowicka’s  equation  (9)  .  The  coefficient  for  /%  can  be  written 


42 

N 


m* 1 

^  xn  , 


_ 


r  -s-  w\ 


\  Y\  -  x  vw  +1  v x  t '  j 
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The  condition  for  equilibrium  for  case  1  is 


V\U 


O  Y 


)tL  -- 
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J 
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Bin 

2Y\-V\ 


_ (,2  8) 


For  the  deflection  of  the  slab  and  the  foundation  to  be  compatible 
equations  (21  and  (27)  must  be  identical.  This  condition  is  fulfilled 
when  the  coefficients  of  all  like  powers  of-£  are  equal.  Writing 


W  AO 
N 


we  have  from  the  coefficients  of 


A 

v\r  r> 


(iVO  +»^V-2V\) 


and  for  the  general  1  erm 


where 


.Yi  - (?<\) 

2.0. 


£ 


v\  =  o 


(x>nn 


-  _X_ 
x.  K 


5 


-  .VJ  — QSp) 

2  wa 


By  splitting 

the  terms  in 

aqua' ions  (29) 

and  (30) 

into  partial 

fractions  and 

substituting 

Bx^ 

J[oW 

c\ 

/ 

Bin 

•  _L 

4ov  VN_ 

V\s  0 

2  r\  vi 

2  m 

'  2.VY1Q 

16. 


we  have  the  following  homogenous  equations  along  with  (28)  for  the 
solution  of  the  b's. 


e> 


a  o 


|-^n 


-V 


yA 


Y\--0 


V\) 


^  0 


- Qm*) 


Eh, 


V\  :  0 


2  YY\  -  2.v\'V 


(2^±)'.  —  o  - f\b) 

^  k  ^2.  ws — \)  \ 


VY\  <:  Z 

Equations  (31)  c?n  he  solved  in  terms  of  B0  and  the  final  expression 
for  the  coefficients  can  be  obtained  in  terras  of  the  load  and  the  half 
width  of  the  slab  from  (28).  The  value  of  the  resulting  coefficients 
substituted  in  the  expression 


^  -  z.  M-C 


gives  the  required  formula  for  the  determination  of  the  contact  pressure 
of  a  concentrated  line  load  at  the  edge  of  the  slab  when  the  ’  5 

are  known. 

Case  2 


Overturning  Moments  at  the  Edge  of  the  Slab 


z 


F  \  0 
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For  the  case  of  the  overturning  moments  shown  in  figure  4  the 
boundry  conditions  for  the  upper  surface  of  the  foundation  are 

X  >  <x  -  o 

/  <  a  --  ^00  . 


For  the  slab  we  have 

Q 


cTqo  -  -  W\. 
d  y-2-  D 


*  =  O 


to  =  uJ0 


The  slab  equation  for  case  2  is 


)  - 
d  ** 


O 


d 


Proceeding  in  a  manner  similar  to  that  for  the  concentrated  loads  it 
results  that  the  deflection  of  the  slab  can  be  expressed  as 


The  deflection  of  the  foundation  is  given  by  equation  (27).  Comparing 
coefficients  of  like  powers  of  we  obtain  from  (27)  and  (32)  for  \flaj 


W  3  0 


fern  h  _3L 
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The  condition  of  equilibrium  requires  that 


I 

V\1  0 


&  X  v-\__  -  O 

2  Y\-*\ 


-  *  -  1^) 


Breaking  the  terms  in  equations  (32)  into  partial  portions  and 
combining  with  (33)  the  B*s  in  terms  of  Bn  can  be  determined  from 
(33)  and  the  following  homogenous  equations 


v\  *  O 
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2m  -xy\-\ 
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(2.  W\-  V)  '• 
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where  vv\  2L  2  The  final  value  of  the  coefficients  in  terms 
of  W,  a  and  W  can  be  determined  by  substitution  in  (35) 
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The  resulting  expression  for  the  coefficients  placed  in  the  equation 


=  \ 


"  X 


2n 


will  give  the  required  formula  for  the  determination  of  the  contact 

pressure. 
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Case  3 
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This  case  has  been  treated  by  Borowicka  and  the  resold  ing 


equations  for  +  he  solution  of  the  B's  are 
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where  ''As  t  By  breaking  the  summation  expressions  in  (37)  into 
partial  |Vqct ions  and  combining  with  the  condition  of  equilibrium  (36) 
we  have  the  following  homogenous  equations  (38)  and  (39)  for  the 


determination  of  the  B's  in  terms  of  B0: 
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- ^ - - 
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(i.WV-0  \ 


The  B*s  are  finally  evaluated  in  terms  of  the  load  ox  from  the 
equation  of  equilibrium  (36).  Substitution  of  the  resulting  coeffi¬ 
cients  in  the  equation 


will  establish  the  expression  necessary  to  determine  the  contact 
pressure  distribution. 

The  labour  require  to  determine  the  coefficients  for  cases  1 
and  2  can  be  reduced  by  realizing  that  the  group  of  equations  given  by 


m  z.  \ 


is  common  to  both  solutions.  A  group  of  these  er nations  for  any  given 
K  value  could  be  solved  in  terms  of  30  and  ^2)  then  by  substitution  \v% 
the  two  remaining  equations,  a  solution  of  the  unknown  coefficients  in  ' 
terms  of  the  load  -results.  The  two  remaining  equations  for  case  1  are 
(31)  (a)  and  (23),  and  for  case  2  (35)  and  (33). 
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Discussion  of  Borowickals  Solution 


21. 


In  attempts  +j0  determine  the  solution  for  the  contact  pressure  of 
symmetrically  placer!  concentrated  or  discontinuous  uniform  loads  it  was 
found  that  sufficient  equations  could  not  bo  obtained  to  compare  all  like 
powers  of  rj  a  It  was  found  that  the  equation  for  the  deflection  of 
the  slab  contained  even  and  odd  powers  of  x  whereas  that  for  the  foundation 
had  only  terms  in  even  powers  of  x  and  a  term  in  log  >/q  * 


To  indicate  the  reason  for  this  discrepancy,  the  integral  involving 
the  deflection  of  the  foundation  is  first  considered.  The  deflection  oPQt) 

is  expressed  as 


=  \  \  \ocj  \  b  -  >'  -t  locj  V  ^  +  *■)  -  Z.  lOCj  i  |  ,  ^ 


expanding  we  have 
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Defining  the  integral 

\  riji 
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we  have  integrating  1-^  by  parts 
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I,  =  "  l  b°)  -  F(f))  'OCJ  t-  .  -V  Fl\)  '  Ftf)y  \  (,V-  *) 


F^l  -  Fw  a  F 

V*  1 


For  the  seond  integral  we  have,  integrating  by  parts 

1,-  \ 

'  n 


1,  F  0)  \ocj  t^  +  x)  -  Ftp)  b<j  x  —  S  F(Sp  ci^j 
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Expanding  Fl^)  in  the  last  term  of  (41)  in  a  Taylor  series  about  (»x) , 
we  have 

F^)  =  Fl-b)  t  Ft:*)  (^t*)  +  Fc-Vj  +  *)  +  -  -  -  - 


Substituting  the  expanded  function  back  in  (41)  we  can  write  that 


Flo  \o c^(l-+Tt) 


which  upon  integration  gives 
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Collecting  terms  in  1^  and  we  have 
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When  dealing  with  loads  that  are  continuous  over  the  full  width  of  the 
slab  and  are  symmetrically  placed,  the  contact  pressure  can  be 

defined  in  even  powers  of  for  x  ^othat  is 


Then  r^)  which  was  defined  as 


can  be  written  in  series  form 


x  ?  o  . 
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F(*i  -  Co  v  c^*.1  't  +  -  -  - 
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In  the  limit  as  we  have  that  the 


\  VW\  o 


Now  considering  equation  (42)  for  w(x)  we  see  that  the  log  x  term 
disappears  and  that  precisely  the  same  terms  which  appeared  in  Boro- 
wicka's  solution  remain.  This  argument  then  establishes  the  validity 
of  the  use  of  a  Maclaurin  series  with  even  powers  to  represent  the  contact 
pressure  distribution  for  continuous  symmetrical  loadings. 

Let  us  now  consider  the  case  of  the  symmetrical  but  discontinuous 
loading.  For  example  the  concentrated  load  at  the  centre  of  the  slab. 

The  contact  pressure  is  now  defined  by  both  even  and  odd  powers  of  x  as 


X  7  o 


We  defined 


or  in  series  form  we  can  write 


Ffr)  -  C*  +  M  Cx  ^  -t 
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In  the  limit  as  >-*■<?  we  have  that 
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An  example  of  a  function  that  has  such  a  discontinuity  is  the 
shear  under  a  concentrated  load.  The  log  x  term  in  equation  (42)  no 
longer  cancels.  This  implies  that  there  is  some  form  of  singularity 
at  the  origin.  The  assumption  then  that  the  contact  pressure  p(x)  can 
"be  expressed  as  a  power  series  containing  even  and  odd  positive  powers 
of  x  is  not  valid.  Therefore  for  discontinuous  symmetrical  loads  a 
different  series  expansion  must  be  considered  or  a  different  method  of 
analysis  investigated. 

It  should  be  noted  that  substitution  of  \-y-  or  y.'--  \Jcf'  in  (42) 
would  not  remove  the  singularity  that  occurs. 

Inverse  Problem 

A  consideration  of  what  is  termed  the  inverse  problem  leads  to  a 
determination  of  what  superimposed  load  should  be  placed  on  the  slab  to 
produce  a  prescribed  deflection  curve  or  a  prescribed  contact  pressure 
distribution.  Or  fundamentally  what  physical  interpretation  can  be 
attached  to  an  assumed  expression  for  the  contact  pressure  or  the  deflection 
curve. 

Let,  us  firs+  consider  the  case  where  the  contact  pressure  is 
assumed  and  is  of  the  form 

=  •*.  1^-j5  .  ---- 

where  s  is  a  positive  integer.  The  deflection  of  the  foundation  can  be 
commuted  from  Bouessesq’s  equation  (12),  and  proceeding  as  under  case  1 

\  \ ocj  i  \-A  -  \ 


we  have  that 
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and  upon  substitution  of  (43)  and  completion  of  the  integration  the 
expression  for  the  deflection  reduces  to 


UJ0  -VJO 


$1± 

IW 


The  terras  in  the  firs4-  square  bracket  of  expression  (44)  apply 
only  when  s  is  an  even  integer,  and  those  in  the  second  bracket  apply 
only  when  s  is  odd.  The  equation  for  the  determination  of  the  super¬ 
imposed  load  q(x)  is  obtained  by  first  differentiating  (44)  with  respect 
to  x  and  then  substituting  in  the  slab  equation 


D  d4 
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As  a  first  example  consider  the  case  when  ip(jq  is  a  constant  « C 
that  is  ■&  =  O  We  have  then  from  (44) 
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Differentiating  (45)  and  summing  the  series  we  have 
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Substitution  of  (46)  in  the  slab  equations  gives  the  expression 
for  the  superimposed  load  on  the  slab  as  follows 


28. 


The  distribution  of  the  superimposed  load  when  the  foundation 
reaction  is  assumed  constant  is  shown  in  figure  6. 


As  a  second  example  consider  the  problem  when  the  contact  pressure  is 


as  shown  in  figure  7. 
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From  equation  (44)  we  have  for  s =  \ 
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By  summing  the  last  two  series  expressions,  equation  (48)  can  be 
written  in  simrler  form  as 


o  **  'wO  — 
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where  aa  -  */  a  Taking  the  fourth  derivate  of  (49)  with  respect  to  aa. 

we  have 


Substituting  in  the  slab  equation  we  have  for  the  required 
superimposed  load 


where  Cv= 


-  sjgt  fi 


I  j  I 


30. 


The  superimposed  load  for  an  assumed  linear  variation  of  the  contact 

pressure  approaches  when  o  and-*®  when  \ 

<x 

A  plot  of  equation  (50)  is  shown  in  figure  3.  This  solution  is  of 

practical  interest  when  combined  with  others  to  give  an  approximation 
to  more  practical  problems  such  as  the  concentrated  load. 
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By  applying  the  principle  of  superposition  to  the  two  foregoing  examples 
we  obtain  such  an  approximation  for  the  concentrated  load.  The  load  on 
the  slab  is  given  as 

—  M 


and  the  graph  of  the  function  is  shown  in  figure  9. 
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Let  us  now  consider  problems  for  which  the  deflection  of  the  dab 
or  foundation  is  some  prescribed  function.  Considering  the  deflection 
of  the  foundation 


^  =  Jk_ 

TCtd 


-\ 


which  arises  from  Boussinesq's  equation  (12),  we  have  on  differentiating 
with  respect  to  x  that 

-  =  Jx  bL  y  d  . 
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Recalling  that 
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we  have  for  the  above  integral 
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as  £  -*  o  we  have  that 
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where  the  symbol  f'  stands  for  the  principal  value  of  the  integral, 

and  c_ "  — - - 


g 

The  solution  of  (52)  is  obtained  by  first  transforming  the  variables 
by  letting 
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We  then  have  for  (52) 
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Expanding  P@)  iV\  ©  in  its  Fourier  expansion,  we  have 
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Similarly  the  function  G  (.<£>)  is  thought  to  be  developed  in  the  interval 
(o,^)  in  the  following  form 

sw \  <j>  G  ~  ^  if\  yy\  j)  -  — ^  s) 


Substituting  (54)  and  (55)  in  (53)  we  obtain 
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that  is  to  m  =  and  P\0  is  arbitary.  This  then  permits  a  determin¬ 
ation  of  the  function  H.©)  with  unknown  coefficients  to  be  made 

in  terms  of  the  known  coefficients  to), 

Let 
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Developing  as  a  sine  series,  we  write 
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The  integral  on  the  r  V\.^,  is  zero  when  m\uand  is 

'  z 

when  m-.o  Therefore  it  follows  ’:hat  the  coefficients  f\  ,A  can  be 

expressed  as 

i*\A  =  _2_  \  (4  [4>]  tiv\$  Vm  r\<j)  \(j)  -  - -($1) 

' 0 

The  unknown  coefficients  thus  determined  in  terms  of  the  known 

coefficients  of  the  function  C*f>)  give  a  solution  of  the  reaction 
function  P(p>) 

As  a  first  illustration  of  this  method  consider  the  case  when  the 
deflection  of  the  foundation  is  constant  that  is 

uo  =  C0  <  0 

From  the  foregoing  work  we  have 
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Equation  (59)  agrees  with  that  given  by  Timoshenko  reference  6  page  93, 
for  the  pressure  distribution  under  a  rigid  strip  which  has  undergone 
a  constant  deflection. 

As  a  second  example  consider  the  -case  when  the  deflection  is 
assumed  parabolic 
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We  have 
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and  from  (53)  we  write  that 
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Equating  coefficients  we  find  that 
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We  can  then  write  for  the  reaction  function 

Y  C  coi©)  *\v\©  ■  =  cosi©  ~v 

where  AQ  is  an  arbitary  constant.  Transforming  to  the  original  variable 
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Taking  as  finite  at  ^  -  )  we  have  that  f\0-  C. x  and  the  resulting 
expression  for  the  contact  oressure  is 
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By  substituting  (60)  and  (61)  in  the  slab  equation  we  have  that 
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and  therefore 


37. 


That  is  the  superimposed  load  and  the  contact  pressure  are  equal 
and  opposite  and  form  respectively  the  upper  and  lower  halves  of  an  ellipse 
as  shown  in  figure  10. 


the  problem  when  the  deflection  of  the  foundation  is  linear  and  symmetrical 
about  the  centreline,  that  is 
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The  "boundary  conditions  for  the  plate  are 
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It  has  been  shown  +hat  the  slope  can  be  expressed  as 
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This  expression  can  be  transformed  into 
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By  expressing  P(0)  in  the  series  form 
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we  have  shown  by  equation  (57)  that  the  unknown  coefficients  Kn 
can  be  determined  when  V\($)  is  known. 
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Therefore 


From  (57)  we  find 
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Substituting  in  (63)  we  have 


p\  ^  a  Cm  ;.  |  c 0^^A\\)<f>  '  CO t>(,A-0^  ^  -_£}  eo^n  +  0  0  -  co^Cn-v)  (6|d^. 

1 


Integrating, the  coefficients  are  given  by 
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For  y\=  \  we  have 
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which  shows  that 


K*  o 


From  the  trigonometeric  functions  we  have 
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therefore  (64)  becomes 
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(65)  gives  that  all  the  coefficients 
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The  reaction  function  can  then  be  expressed  as 
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To  evaluate  the  series  expression  of  (67)  we  consider  the  function 
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and  select  only  the  real  part  of  the  answer.  In  (68)  we  have  that 
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Let 
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and  consider  each  function  separately.  T?e  have  then 
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To  sum  the  series  given  by  (69)  we  multiply  both  sides  by 
and  differentiate. 
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The  resulting  geometeric  series  can  be  summed  as  follows! 
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To  integrate  the  above  expression  the  l-  Vv*-  is  split  into  partial 
fractions  and  we  find 


Integrating  we  have 
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similar  manner  we  find  for 
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the  resulting  expression  becomes 
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Adding  (71)  and  (72)  7*e  have 
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Recalling  that 
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we  have  on  substitution  that 
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The  art-Van  2  can  be  written  as 
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Therefore  (73)  becomes 
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In  complex  variable  the  principal  value  of  the  logarithmic  function  is 
given  as 
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and  thereby  the  function  | (zj  can  be  written  as 
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Selecting  only  the  real  part  of  (74)  we  have  for  the  reaction  function 
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In  terms  of  the  variable  x  the  foundation  reaction  for  linear 
displacement  becomes 
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The  determination  of  the  constant  (\0  will  depend  on  the  contact 
pressure  at  the  edges  of  the  slab.  When  the  contact  pressure  is  assumed 

i 

finite  at  the  edges,  we  have  from  (7 6)  that  l\o  must  be  zero  and  therefore 
the  foundation  reaction  for  this  case  becomes 
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From  the  slab  eq^iation  we  have  that  =  IpW 

This  solution  gives  an  approximation  to  the  concentrated  load  as  shown 


load  are  equal  and  opposite  and  independent  of  the  physical  properties  of 
the  slab  and  the  foundation  for  all  deflection  curves  of  the  cubic  order 
or  less.  For  deflection  curves  higher  than  the  cubic  the  superimposed  load 
is  dependent  on  the  physical  properties  of  the  slab  and  foundation  and  is 
given  by  the  following  expression 
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Conclusions 

It  has  been  shown  that  Boro /icka*s  method  of  analysis  is  valid  only 
for  a  special  class  of  problems.  The  nature  of  the  problem  must  be  such 
that  the  deflection  function  and  all  its  derivatives  be  continuous.  This 
means  that  the  slope,  moment,  shear,  and  load  functions  cannot  have  a 
singularity  at  any  point  over  the  width  of  the  slab.  It  can  be  readily 
seen  that  such  stringent  conditions  limit  the  applications  of  the  method 
in  actual  practice.  To  solve  the  problem  which  has  a  discontinuous  loading, 
for  example  the  concentrated  load,  or  discontinuous  uniform  load  on  the  slab 
other  scries  expressions  for  the  contact  pressure  and  deflection  curve 
would  have  to  be  assumed;  or  other  methods  of  analysis  investigated.  An 
approximate  solution  for  concentrated  loads  is  obtained  as  suggested  under 
the  part  "Inverse  Problem".  This  particular  method  could  be  extended  and 
a  variety  of  physical  problems  could  be  built  up  by  applying  the  principle 
of  superposition. 

The  solution  of  equation  (52)  for  known  deflection  curves  would  give 
the  contact  pressure  distribution  beneath  actual  structures  such  as  locks. 

A  number  of  these  studies  made  for  various  types  of  soils  and  size  of 
structures  could  be  useful  in  a  preliminary  investigation  of  lock  design. 

For  experimental  work  done  in  a  laboratory,  equation  (52)  would  give  a 
means  of  checking  computed  and  observed  contact  pressures  beneath  elastically 
supported  strips.  The  computed  contact  pressure  being  obtained  from  the 
observed  deflection  curve. 
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In  general  the  solution  of  the  slab  resting  on  elastic  foundation 
would  require  the  solution  of  the  following  integral  equation 
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Expression  (77)  is  obtained  by  the  coalescence  of  the  slab  equation 
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and  the  expression  for  the  deflection  of  the  foundation 
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Another  method  of  obtaining  an  exact  solution  to  the  problem  of 
the  discontinuous  load  would  be  from  a  consideration  of  (52).  It  has 
been  shown  that 

14a  =  A  {(p)  ■-  1  L  Pie)  VvftQ  de  -  -  -  h&J 
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We  kno^  that  the  solution  of  (78)  can  he  obtained  when  is  expressed 

in  the  series  form 
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and  the  contact  pressure  is  given  as 
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A  solution  of  the  foundation  reaction  can  then  be  obtained  from  a 

c  nsideration  of  the  slab  equation  in  which  each  term  is  expressed  by  its 

corresponding  series  expansion.  That  is,  the  superimposed  load  on  the 

slab  is  given  by  a  Fourier  series  with  known  coefficients  C-™  and  the 

fourth  derivative  of  the  deflection  and  the  contact  pressure  are  similarly 

expressed  in  series  form  with  unknown  coefficients  (\ws  •  From  the 

resulting  expression  the  equating  of  coefficients  permits  a  determination 

of  the  unknown  coefficients  to  be  made  in  terms  of  the  known  coefficients 

and  subsequently  a  solution  to  the  problem  is  obtained. 

The  solution  of  the  continuously  loaded  strip  as  given  under  cases  1, 

2,  3,  for  which  PW  is  assumed  to  be  a  power  series,  can  be  determined  with 

lass  effort  from  a  consideration,  of  the  integral 
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than  tha  solution  obtained  from  the  equation  given  by  Borowicka  as 


Equation  (79)  can  be  placed  into  a  more  useful  form  by  writing 
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This  expression  reduces  to 
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The  solution  of  the  contact  pressure  for  the  elastic  strip  on  an 
elastic  foundation  as  given  under  cases  1,  2  and  3  can  be  applied  also 
to  indeterminate  structures  such  as  culverts.  The  moments  at  the 
corners  of  the  culvert  could  be  determined  by  using  the  method  of 


moment  distribution 
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